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TECHNICAL PUBLICATION 


MAGNETIC FLUX COMPRESSION REACTOR CONCEPTS 
FOR SPACECRAFT PROPULSION AND POWER 

1. INTRODUCTION 


Driven by the desire for fast, efficient interplanetary transport, NASA has set a number of ambi- 
tious mission goals for future space flight. These goals include accomplishing missions to Mars within 
months rather than years and improving propulsion capability by a factor of 10 within 1 5 yr. These criteria 
can be roughly translated into a specific power goal of 10 kW/kg for the onboard powerplant and a specific 
impulse (/ ) goal of 10,000 sec in combination with a thrust-to-mass ratio (acceleration) in excess of 
10 -3 g for the spacecraft propulsion system. The principal technical challenge is the development of pro- 
pulsion systems exhibiting rapid acceleration and very high / attributes. 

Unfortunately, most highly efficient propulsion systems which are within the capabilities of present 
day technologies are either very heavy or yield very low impulse, such that the acceleration time to final 
velocity is too long to be of lasting interest. One exception, the nuclear thermal thruster, could achieve the 
desired acceleration but would require inordinately large propellant-to-payload mass ratios to reach the 
range of desired final velocities. 

A particularly promising alternative approach, among several competing concepts that are beyond 
our modem technical capabilities, is a low-yield thermonuclear pulse thruster. 1-14 In this scheme, neutron- 
lean microfusion detonations form an expanding diamagnetic plasma cloud that compresses the magnetic 
flux within a semienclosed reactor structure. The plasma expansion is ultimately reversed by increasing 
magnetic pressure, and the detonation products are then collimated and expelled by a magnetic nozzle. The 
charged-particle expansion velocity in these detonations can be on the order of 10 6 -10 7 m/sec and if 
effectively collimated by a magnetic nozzle, can yield the / and acceleration levels needed for practical 
interplanetary flight. 

Methods for inductively extracting electrical power from the compressed magnetic field can 
also be envisioned. This is an integral component of the scheme since the energy needed to power the 
standoff drivers for the subsequent detonation is extremely high. Furthermore, an indirect mode of opera- 
tion is conceivable in which a detonation-driven electric power reactor is used to drive high-power electric 
thrusters. 

Magnetic flux compression concepts merit reassessment for spacecraft applications, based on modem 
developments in target and reactor design principles. For example, these schemes could capitalize on 
recent advances in target technologies including inertial confinement fusion, 15 staged microfusion 


schemes , 4 - 14 ’ 16 ’ 17 magnetized target fusion , 18 ’ 19 and fast ignition antimatter-initiated nuclear detonations. 20- 
22 It is also conceivable that chemical detonations utilizing advanced high-energy density matter (HEDM) 
charges could attain limited success if the resulting plasma clouds demonstrate the requisite diamagnetic 
properties. 

Moreover, inductive storage power reactors based on field compression processes have been pro- 
posed as a prime pulse power source for high current standoff drivers and as a direct conversion topping 
cycle for pulsed fusion systems . 23-28 In comparison with capacitor banks, inductive energy storage 
devices are generally more compact and less costly and avoid the need for initial opening switches in the 
pulse-forming line. The possible utilization of high-temperature superconductors as a flux compression 
reactor surface appears to offer attractive advantages as well. 

The objective of this work is to assess the system-level performance and engineering challenges of 
field compression reactors as applicable to interplanetary spacecraft. In particular, special attention is 
devoted to devices which compress the field between a detonation plasma armature and a superconducting 
stafor. 
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2. FIELD COMPRESSION REACTOR CONCEPTS 


Magnetic flux compression reactors are based on a multistep energy conversion process as follows: 
chemical/nuclear — > kinetic -> electrical -» kinetic. The initial detonation charge is first transformed into 
the kinetic energy of a moving conductive armature. This kinetic energy is then transformed into electrical 
energy as the armature compresses an initial seed magnetic flux against a highly conductive stator surface. 
More specifically, a fraction of the energy is temporarily stored in a rapidly increasing magnetic field while 
a portion is extracted from the induced stator current as useful electrical power. If the kinetic energy of the 
armature is fully expended without destroying the device, then the compressed magnetic field can reverse 
the armature motion and return the stored electromagnetic energy as kinetic energy of the armature. The 
reaccelerated armature can then be expelled from the device for direct thrust production. Of course, all of 
these processes are subjected to conversion inefficiencies associated with nonideal losses. 

In reactors of this type, magnetic flux is trapped between the fast-moving armature and the station- 
ary stator because of circulating currents that are induced within the field compression surfaces. In fact, the 
magnetic flux can be trapped only if the armature and stator have electrical conductivities that are high 
enough to resist magnetic field penetration during the field compression process. This criterion naturally 
leads to consideration of highly conductive detonation sources (i.e., fusion plasma) and highly conductive 
stator materials (i.e., superconductors). 

From the viewpoint of global energy conservation, these reactors obey the relation: 

W jet + W ekc = W k -W a -W s = T] k W k , (1) 

where W et is the propulsive jet energy of the exhaust, W e j ec is the useful extracted electrical energy, is 
the initial kinetic energy of the detonation charge, W a is the energy dissipated in the armature, W s is the 
energy dissipated in the stator circuit, and is a global energy conversion efficiency. 

Ideally, the energy conversion efficiency can exceed 80 percent, but this value could be reduced 
substantially by real plasma processes such as flute instabilities, electron Joule heating effects, and field- 
aligned ion flow due to the ambipolar potential. Furthermore, the magnetic field configuration and reactor 
chamber design plays a significant role in terms of magnetic diffusion losses and magnetic flux compres- 
sion efficiency. 

For a pure electric power reactor, one is interested in achieving the highest possible W elec with little 
concern for propulsive performance. For a propulsive reactor, one is interested in achieving the highest 
possible W et while also extracting the minimum necessary electrical energy that is needed to ignite the 
subsequent* detonation. Two reactor design concepts suitable for spacecraft applications are discussed in 
section 2.1. 
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2.1 Propulsion Reactor 


The typical propulsion reactor concept is based on a concave surface around which magnet coils 
are situated to form a magnetic mirror. In addition, because most pulse fusion schemes require a large 
amount of driver energy, these reactors almost always include a generator coil for electrical power genera- 
tion. The Daedalus-class magnetic reaction chamber exemplifies this highly integrated approach. 6 * 7 In this 
technical paper, a slight variation on the Daedalus theme in which a bulk-processed, high-temperature 
superconductor acts as the reaction chamber wall is considered. This concept is illustrated in figure 1. 



Figure 1 . Magnetic flux compression reactor concept for integrated 
space propulsion and power. 

In this concept, the plasma armature expands against a rearward-diverging magnetic field provided 
by a set of superconducting coils. The magnetic flux is trapped and compressed between the plasma and a 
high-temperature superconductor (HTSC) reaction chamber wall (i.e., stator). At some point, the expan- 
sion process is reversed and the increased magnetic pressure acting on the plasma forces it out the rear of 
the reaction chamber. The compressed magnetic field also produces a magnetic pressure on the reaction 
chamber wall, yielding a forward-thrust component. Electrical power is extracted from the flux compres- 
sion process using an inductive pickup coil located at the exit of the reaction chamber. 

By introducing an HTSC stator, it should be possible to reduce flux diffusion losses, reduce the 
minimum seed field, and reduce the overall size of the reaction chamber. The hoped-for result is a compact, 
low- weight integrated propulsion and power reactor driven by low-yield nuclear fusion detonation pulses 
that can yield acceptable interplanetary trip times while carrying sufficient payload to perform meaningful 
missions. 

A nuclear pulse system based on this scheme can provide performance levels in the required regime 
of both high jet power and high / . Along these lines, it is instructive to compute the theoretical perfor- 
mance limits that can be achieved for various fusion fuels. 
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2.1.1 Energetics 

Starting with the assumption that after the nuclear reactions are quenched by the hydrodynamic 
disassembly, the fusion energy is averaged by collisions among all the particles in the plasma including the 
blowoff mass which did not undergo reaction. After complete thermalization, the energy per ion is the 
same throughout the fireball and the following relationship for energy conservation can be written: 


1 2 /fcAmc 

— m;U; = 

2 N; 


( 2 ) 


where m i is the mass of ion /, u i is the velocity of ion ij b is the bum fraction. Am is the total mass defect for 
the fusion reaction, c is the speed of light, and N t is the number of ion species generated per reaction. 
Neutrons are neglected because they carry no charge. Only the momentum of charge particles can be 
electromagnetically manipulated for direct thrust production. By introducing a parameter for the mass 
fraction that is converted to energy by the reaction (cc= Am/m R ), it is possible to write eq. (2) in the form. 


(3) 

2 JV; 

where m D = m . From momentum conservation considerations, the following relationship is obtained: 

A A l 


"W = \[{ l - fbTL R m i U i + fbL P m i u \ ’ (4) 

where the summations are over the reactants ( R ) and the products ( P ) and the overbar denotes an average 
over all ion species. Equation (3) can now be used to eliminate on the right-hand side of eq. (4), such that 
an expression for the effective ion expansion velocity as a function of the bum fraction. 


u i ~ [O “ + fbH P 4™i] J 2 -. 

where the mean ion mass is defined by 


f b am R c 


N: 


(5) 


"H =^[( i ~fb)I 4R m i + fbLp m i 


The theoretical I sp may now be computed by noting that 


. / m.Uj _ 

I S p = = = -= L = u i 

F mi nil 


( 6 ) 


(7) 
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or when referenced to the Earth’s gravitational field, 



Ui_ 

SO 


( 8 ) 


where g 0 = 9.81 m/s 2 . 

The nuclear fusion reactions of practical interest involve isotopes of hydrogen such as deuterium 
(D) and tritium (T) and isotopes of helium (He). These reactions are given in table 1 where the subscripts 
denote the nuclear charge (i.e., number of protons) and the superscripts denote the atomic mass number 
(i.e., total number of nucleons). The proton and neutron branches of the D-D reaction have similar valued 
cross sections such that each branch consumes roughly 50 percent of the fuel. The energy yield parameters 
for the various fusion reactions are given in table 2. It is important to note that the D-T and D-D reactions 
have larger fusion cross sections than the D-He reaction and are easier to ignite; however, they release 
approximately 80 percent and 35 percent, respectively, of the reaction energy in the form of neutrons 
which are useless for flux compression and wasteful for space propulsion applications, since they cannot 
be manipulated by electromagnetic means. 


Table 1. Nuclear fusion reactions of major practical interest. 


D-D 

1 D 2 + 1 D 2 -» ,13(1 .01 MeV) + ^(3.02 MeV) 

[W= 4.03 MeV] 


-» 2 He 3 (0.82 MeV) + 0 n' (2.45 MeV) 

[W=3.27 MeV] 

D-T 

,D 2 + ,T 3 -> 2 He 4 (3.5 MeV) + on 1 (14.1 MeV) 

[W= 17.6 MeV] 

D-He 

1 D 2 + 2 He 3 -> 2 He 4 (3.6MeV) + 1 p 1 (14.7MeV) 

[147=18.3 MeV] 


Table 2. Energy yield for major fusion fuels. 


Fuel 

Reaction 

Products 

w(J/kg) 

a 

D-D 

iT 3 , iP 1 

9.9 xIO 13 

i.i xi<r 3 

D-D 

2 He 3 , on 1 

8.1 xIO 13 

CO 

0 

X 

1 

D-T 

2 He 4 , on 1 

3.38x10 14 

3.75 x10“ 3 

D-He 

2He 4 , 1 p 1 

3.52 xIO 14 

3.91x10- 3 


Practical considerations compel the use of D-He as a fusion fuel for space propulsion. Although the 
required ignition energy is substantially higher for D-He, the reaction products will consist predominately 
of charged particles, which can be electromagnetically manipulated. Of course, there will always be some 
neutron production due to D-D reactions of the primary fuel and D-T reactions between the primary D and 
secondary T. However, calculations indicate that the number of neutron reactions can be <5 percent of the 
D-He reactions. 6 
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Using the parameters summarized in tables 1 and 2, it is possible to compute the effective detona- 
tion expansion velocity and theoretical I sp for the D-He fusion reaction of interest using eqs. (5) and (8). 
The results of these calculations are shown in figure 2 as a function of the bum fraction. Detailed inertial 
confinement fusion calculations indicate that it should be possible to achieve a bum fraction of at least 10 
percent, yielding an effective expansion velocity approaching 10 7 m/sec. The I sp in this range will be on 
the order of 10 6 sec. Note that these calculations represent only a theoretical upper limit. 



Burnup Fraction, ft 



Figure 2. Computed effective expansion velocity and I sp for a D-He fusion detonation. 
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2.1.2 Thrust Performance 


The average thrust (7) produced by a nuclear pulse engine is governed by the pulse rate (f p ), the 
pellet mass ( m p ), the magnetic nozzle efficiency (r| j), and the effective expansion velocity as defined by the 
following relationship: 


T=rijf p m p u i . V) 

The average thrust of a propulsive reactor has been calculated assuming a 200-MJ yield fusion 
detonation of D-He, a pulse rate of 100 Hz, and a nozzle efficiency of 65 percent. The results are shown as 
a function of bum fraction in figure 3. When assuming a constant energy yield, the average thrust decreases 
with increasing bum fraction, since the pellet mass is decreasing with increasing bum efficiency. By the 
same token, the I sp will increase with increasing burn efficiency. At a practical bum fraction of 10 percent, 
the average thrust is ~3 kN and the / is approximately 8x 10 5 sec for the assumed parameter values. 



Figure 3. Average reactor thrust versus bumup fraction assuming 200-MJ 
fusion detonations of D-He at a pulse rate of 100 Hz with a 
nozzle efficiency of 65 percent. 
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2.2 Power Reactor 

A schematic of a simple radial-mode power reactor with inductively coupled circuits for power 
extraction is shown in figure 4. Although other geometries are conceivable (e.g., spherical), the radial 
reactor is more amenable to power production with neutralized thrust. In the radial configuration, the 
generator consists of three nested shells. The innermost shell is the pickup (stator) coil in which electrical 
power is inductively extracted from the flux compression process. It is anticipated that the pickup coil will 
be fabricated from HTSC material. The pickup coil is surrounded by a conductive shielding shell which 
may or may not be fabricated from an HTSC material in bulk processed form. The purpose of the shielding 
shell is to minimize induced current transients in the outer magnet coil whose function is to produce the 
initial seed flux in the generator. 


Shielding Magnet 




i 


Excitation Magnet Pickup Coil 


Circuit Circuit 





Figure 4. Schematic of radial-mode magnetic flux compression power reactor 
with inductively coupled power extraction circuit. 


If a nuclear detonation is used as the plasma source, a target pellet would be fired into the reaction 
chamber and ignited by a high-energy driver (i.e., laser beam or particle beam). If a chemical detonation is 
used as the plasma source, the explosive component could either be fired into the reaction chamber as a 
target pellet or directly inserted as an expendable cartridge. 

In practice, all of the individual circuit elements will interact inductively. The coupling of these 
circuits is defined by the mutual inductance between the various elements. In a good engineering design, 
the mutual inductance between the excitation magnetic coil and the other circuit elements in the system 
will be negligible. For good energy conversion efficiency, the goal is to cause an effective reduction in 
generator self-inductance through an increase in the mutual inductance between the expanding armature 
and the pickup coil. In this way, it is not necessary to deform the circuits. 
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As it is envisioned, the power reactor would embody superconducting magnet energy storage (SMES) 
technology as a means of inductively storing the extracted electrical power. This approach, illustrated in 
figure 5, is much more compact, lower in weight, and lower in cost than capacitive energy storage tech- 
nologies. In this system, the load coil would be a normal metal conductor while the storage coil would be 
a superconductor. The SMES coil would be open during pulse generation and immediately closed at peak 
current. An operating reactor must be expected to supply its own operational power requirements, and this 
system would be able to supply the on-demand pulse power needed for ignition of the subsequent detona- 
tion. If the electrical energy is to be used in a normal power grid, several storage coils could be incorpo- 
rated into the design such that continuous direct current power can be conditioned and bussed throughout 
the spacecraft. 


Armature 



Magnet Pickup Load 

Coil Coil Coil SMES 


Figure 5. Schematic of SMES circuit coupled to a flux compression power reactor. 


2.2.1 Simplified Circuit Analysis 

The flux compression reactor circuit is illustrated schematically in figure 6. For simplicity, the 
generator is depicted as a single-turn induction coil connected to a single-turn load coil. The induction coil 
serves as a stator while the expanding plasma serves as an armature. The circuit diagram for the generator 
is shown as a time- varying generator inductance (LJ connected to a fixed load inductance (L l ) with a 
resistive loss component ( R ). 


10 


l"l ' I I 111 HI Mllll'linWXIBUlIMM II 





Ll 


Figure 6. Schematic of flux compression reactor electrical circuit. 


The principle of operation is as follows: (1) Seed flux injection, (2) ignition of centrally located 
charge, (3) detonation-driven expansion of the plasma armature, (4) magnetic flux trapping between the 
plasma armature and the stator coil with magnetic flux diffusion into both the stator and armature, (5) rapid 
field compression as the kinetic energy of the plasma is transformed into magnetic pressure, and (6) forced 
reduction in generator inductance yielding a fast rise time current pulse. 

The circuit equation for the generator may be written in differential form as 


—{U) + Ri = 0 , 
dr 

where L = L + L, and i is the instantaneous current. Solving for the current, the following 
is obtained: 


( 10 ) 


i = ^S.exp{-J'(ff/i)di| (11) 

and find that the current varies inversely with the generator inductance. It is convenient to rewrite eq. (1 1) 
in terms of a flux coefficient A: 

A = = -pr- = exp j-J* (A / L)dt \ < 1 , (12) 

00 ^O'O f 0 J 

where (j) is the magnetic flux. Thus, A is a measure of flux compression efficiency and quantifies the 
magnetic flux diffusion losses into the stator and the armature through a negative exponential function. 
These magnetic diffusion processes are illustrated in figure 7. It should be noted that the expression for the 
flux coefficient demonstrates that flux diffusion losses can be effectively expressed in terms of a resistive 
circuit component. 
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Figure 7. Illustration of magnetic diffusion processes in the stator 
and armature during field compression. 


An expression for the instantaneous current is obtained from eq. (12) in the form: 


'0 


(13) 


The initial energy in the load is W 0 = Vi L 0 / 0 2 and the energy delivered to the load during field compression 
is W L = VlLJ 2 . Therefore, by taking the ratio of W L to Wq and substituting for / using eq. (13), the follow- 
ing expression is obtained: 


{ L A L ) 


w 0 UAl 


(14) 


The power delivered to the load P L is defined by 




(15) 


We note, however, that eq. (10) can be written as 


f ->+4). A - 


( 16 ) 


and that W, =V iLA 2 . Thus, the delivered power takes the form 


P L =-2W L {R+L g )[L 


(17) 
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The major point to note in the preceding circuit analysis is the impact of magnetic diffusion losses 
on generator performance. For example, the current multiplication ratio varies in proportion to A (i/t 0 ~ A), 
and the energy multiplication ratio varies in proportion to A 2 (W L /W 0 ~ A 2 ). Therefore, if one hopes to 
achieve good reactor performance, it is essential that diffusion losses be minimized. This requires that both 
the stator and armature possess low magnetic diffusivity properties. 

The fundamental definition of magnetic diffusivity in terms of material properties is 


Dm ~ 



( 18 ) 


where fi is the magnetic permeability and <7 is the electrical conductivity. Thus, it is desirable to utilize 
materials having the highest possible electrical conductivity. 
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3. MAGNETIC DIFFUSION ISSUES 


The impact of magnetic diffusion losses is equally important to the performance of both propulsive 
and power reactors. The issues associated with stator and armature diffusion losses, including suitable 
design strategies, are discussed separately in sections 3.1 and 3.2. 

3.1 Armature Concepts 

The conventional approach in sacrificial magnetic flux compression reactors is to utilize a metal 
compression surface to minimize diffusion losses. In almost all cases, these designs incorporate a metal 
armature which is deformed and accelerated by high explosives. Magnetoimplosive generators, for 
instance, utilize this design principle with great success. 

In principle, a target/cartridge is considered that is configured as a detonation charge surrounded by 
a thin metal liner, as illustrated in figure 8(a). However, the major drawbacks for the explosively driven 
metal liner ap proach a re as follows: (1) The limited armature speed that can be achieved, (2) the likelihood 
of self-destruction which negates the desire for intermittent firing capability, and (3) general incompatibil- 
ity with microfusion explosion schemes. For example, if one desires high-energy output (say 10^ J) over a 
very short time interval (say 10 jj sec) large-scale devices are required. This implies the need for a rapidly 
accelerated armature, a feat that is difficult to accomplish with existing technologies. 




Detonation 

Plasma 


(a) (b) 

Figure 8. Illustration of explosively driven armatures: 

(a) metal-lined cartridge and (b) detonation plasma. 


An obvious solution to this dilemma is to attempt the use of a detonation plasma armature, as 
illustrated in figure 8(b). By eliminating the metal liner and relying on the intrinsic diamagnetic character- 
istics of the plasma cloud, it may be possible to obtain the desired field compression rate in an intermit- 
tently fired device. A comparison of metal versus plasma armature characteristics is shown in table 3. 
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Table 3. Metal versus plasma armature characteristics. 


Parameter 

Metal 

Plasma 

Conductivity 

Very High 

Relatively Low 

Temperature 

Very Low 

Very High 

Velocity 

Very Low 

Very High 


The major advantages associated with the use of a plasma armature rather than a metal armature 
may be summarized as follows: (1) Greater speed (shorter pulse), (2) greater expansion/compression, (3) 
lower impulse delivered to reactor structure, (4) intermittent firing capability, (5) exhaust products can be 
readily utilized for propulsion, and (6) lower cost. These advantages are of special significance to space 
flight where simplicity and robust performance attributes are essential. On the other hand, plasma armatures 
introduce substantial technical risks that must be addressed and overcome through research and develop- 
ment. 


The major technical risks/uncertainties associated with plasma armatures are as follows: (1) Achiev- 
ing sufficiently high electrical conductivity to ensure diamagnetic behavior, (2) electron Joule heating 
effects, (3) field-aligned ion flow due to the ambipolar potential, (4) Rayleigh-Taylor surface instabilities, 
and (5) assurance of armature rebound. 

The diamagnetic properties of the plasma armature can be quantitatively measured by the relative 
value of the magnetic Reynolds number, a parameter which naturally evolves from dimensional analysis of 
the magnetohydrodynamics (MHD) equations in the form 

B’ 

Rm=-^ nd - = ^0 (J 0 u 0 L 0 ’ (19) 

&app 

where 0 subscript denotes an appropriate characteristic value for the problem. The plasma electrical con- 
ductivity is <7 0 , Uq is the plasma velocity, and L 0 is a characteristic length of the reactor (e.g., diameter). In 
a general sense, R m may be thought of as a relative measure of the induced magnetic field in the plasma 
(B jnd ) with respect to the external applied magnetic field (B app ). Because the magnetic Reynolds number is 
inversely related to the magnetic diffusivity D m 


R, . (20) 

D m 

It is clear that for magnetic Reynolds numbers greater than unity, where the induced field is at least 
as strong as the applied field, the plasma will be resistive to magnetic diffusion. That is, the induced field 
associated with eddy currents in the plasma act in an opposing direction to the applied field, and if the 
induced field is high enough, the plasma interface elfectively behaves as a magnetic compression surface. 
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The important point to note is that the magnetic Reynolds number depends primarily on the product 
<7u such that one would desire high electrical conductivity and high expansion velocities in order to achieve 
low-flux diffusion losses. For fusion detonations, both of these parameters are extremely high and R m »1, 
as desired. For conventional chemical detonations, the characteristic values are marginal and R m > 1 can be 
achieved using carefully designed plasma-jet sources. It is believed that detonation plasmas using 
advanced HEDM materials may achieve R m > 1 , although this hope has yet to be demonstrated in the open 
literature. 


3.2 Stator Concepts 

The simplest approach to stator design is to utilize a metal having a high electrical conductivity, 
such as molydenum. In this way, one can satisfy the demands for both low magnetic diffusivity and high 
material strength in a relatively straightforward manner. Alternatively, possible utilization of a type II 
HTSC reaction chamber to confine the magnetic field as it is compressed outwardly by the expanding 
plasma is suggested. A type II HTSC stator, although prone to substantial technical risk, should be able to 
substantially reduce magnetic diffusion losses and significantly improve energy conversion efficiency. 

Superconductors are divided into two types, depending on their characteristic behavior in the pres- 
ence of a magnetic field. Type I superconductors comprised of pure metals tend to repel a penetrating 
magnetic flux due to the Meissner effect. In the Meissner state, the penetrating field is completely repelled 
from the interior of the superconductor. 

For type I superconductors, the material is either perfectly conducting and exists in the Meissner 
state or it undergoes a phase transformation to the normal state. The major obstacle in using type I super- 
conductors in a flux compression application is the low threshold values (i.e., temperature, current density, 
and magnetic field strength) defining critical transition to the normal state. The maximum critical field for 
a type I superconductor is ~0.2 T, which is far too low for practical application in flux compression 
devices. . 

Ceramic-based type II superconductors, on the other hand, are potentially useful for flux compres- 
sion applications due to the high critical threshold values which can be obtained. For example, the critical 
temperature is above liquid nitrogen temperatures over a wide range of current densities, and superconduc- 
tivity can persist under applied magnetic fields exceeding 100 T. 

In a type II superconductor, however, there are two critical field levels. The first critical field (B c j) 
defines the limiting value for maintaining a true Meissner state and is normally very small. When the 
applied field exceeds B c j, the material enters the so-called mixed state where the field penetrates in quan- 
tized amounts of flux. These points of penetration, known as fluxoids, may be envisioned as circulating 
vortices of current. The second critical field level (B c 2 ) defines the transition to the normal state, and it can 
be large enough for practical flux compression application. Flux penetration in both the Meissner state and 
the mixed stat e as we ll as the characteristic variation of the critical fields for a type II HTSC are illustrated 
schematically in figure 9. 
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Figure 9. Illustration of superconductor field penetration in the Meissner state and the 
mixed state. Characteristic variation of the critical fields for a type II HTSC. 


The interaction of the fluxoids with defects in the superconductor alters its conductive properties; 
that is, a fluxoid encompassing or adjacent to a defect has its energy altered and its free motion through the 
superconductor inhibited. This phenomenon, known as flux pinning, causes a field gradient in the super- 
conductor and gives rise to a net current in the material. Since the pinning force is small, fluxoids can be 
broken loose from their pinning centers, resulting in a net creep of the flux through the conductor as a 
function of time. This results in an effective voltage in a type II superconductor. 

If the current density is low and the magnetic field is not intense, flux creep is insignificant and the 
induced voltage and effective resistance will be essentially zero. At very high fields and high current 
densities, fluxoids will migrate rapidly, giving rise to a phenomenon called flux flow. The effective resis- 
tance can be non-negligible in the flux flow case, and breakdown to the normal state can ensue. This can be 
a particularly exacerbating issue when the flux skin depth is small and the induced current is high; because 
under these circumstances, the current density can be very high and may exceed the critical value for the 
HTSC. 


Note that the magnetic diffusivity of a type II superconductor can be several orders of magnitude 
below that of the best metals, even in flux flow mode. For instance, recent studies investigating the mag- 
netic diffusivity characteristics of various superconductor materials under pulsed magnetic fields indicate 
that the diffusivity can range from 10~ 2 to 10 -6 m 2 /sec. 29 ’ 30 This range can be compared with nominal 
values for various metals as indicated in table 4. 


Table 4. Electrical and magnetic diffusivity properties of common metals. 


Parameter 

Dimension 

Copper 

Brass 

Aluminum 

Stainless 

Steel 

Po (20 °C) 


1,7 

6.2 

2.8 

72 

o 0 (20 °C) 

10 6 (ft-m)" 1 

63.3 

15.7 

39.2 

1.38 


m 2 /sec 

1.26 

5.10 

2.04 

58 
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It is clear that for a type II HTSC, flux penetration will occur under an intense pulsed magnetic 
field. Therefore, some of the major questions that need to be answered are as follows: (1) How resistant are 
HTSC materials to magnetic field penetration in the mixed state? (2) How good are the magnetic diffusion 
characteristics in the mixed state? and (3) How does the HTSC material behave under intense applied- 
pulse fields, and will it break down? 

The hope is that the transient magnetic field will exhibit a penetration time longer than the charac- 
teristic pulse time of the armature expansion process. This concept is illustrated schematically in figure 10 
where a section of the HTSC stator wall is shown with an increasing magnetic field H inside a radial-mode 
reactor. The anticipated magnetic field profiles in the stator wall are shown at various instances in time 
during the armature expansion process. Note that there is some characteristic time (At p ) defining the pen- 
etration time of the magnetic field through the stator and some characteristic time defining the expansion 
time of the detonation plasma (r D ). For practical purposes, it is necessary that At p >x D . 


Superconductor 



Figure 10. Illustration of transient magnetic field diffusion 
and penetration through an HTSC stator. 


The HTSC stator concept is entirely hypothetical at this point and there are a whole host of signifi- 
cant research issues that must be addressed prior to practical implementation. The major uncertainties may 
be categorized as follows: (1) Breakdown of HTSC under strong pulsed fields, (2) hysteresis cycling of 
magnetization, (3) Joule/neutron heating of the material, (4) structural integrity under cyclic loading, and 
(5) bulk processed versus wire fabrication. 


j 
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4. FLUX COMPRESSION DYNAMICS 


Accurate analysis of the reactor requires formal solution of the MHD equations; however, it is 
instructive to develop a simplified analytical model for describing the dynamics of flux compression. This 
approach is adopted here, and a simple model for the radial-mode reactor developed (illustrated in fig. 4). 
This analysis includes finite conductivity effects and attempts to quantify flux diffusion losses and arma- 
ture rebound criteria. The modeling approach is described in sections 4.1 through 4.4. 

4.1 Field Amplification 

During the explosive expansion of the plasma armature, some fraction of the trapped magnetic flux 
will diffuse into the armature and the stator. The flux which escapes from the annular containment region 
is lost for further compression and represents an inefficiency in generator performance. The field amplifi- 
cation, including losses, may be expressed in terms of the flux coefficient as defined in eq. (12): 


0 B \ 


, = *( 

’ 2 2\ 
r s ~ r a J 

(21) 

<Po Bq 

{nrj-nrf- 

) Ho 

/ 2 2’ 
^ ~ r i 

) 


or 



where r is the internal radius of the magnetic flux containment stator, r ( « r s is the initial radius of the 
armature, and r a (t) is the time-dependent radius of the plasma armature. 

Differentiation of eq. (22) yields 


1 dH -ni 

H 0 dt 


r a{ r } ~ *? 2 ) 

(t-i ) 2 


,. 2 \ 


+ 


\ r s 


n 


r a J 


dX_ 

dt 


where u a = dr a / dt. Using eq. (22), the following is obtained: 


(23) 
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( 24 ) 


±™ = 2K JL i_^ + 

Hq * ffo% 2 -r„ 2 ) 


H 1 

dX 

Hq A 

dt 


or 


1 dH 2r a 

1 

d A 

H dt ~r?-r} Ua 

A 

dt 


(25) 


This result demonstrates that the maximum field amplification can occur during armature expansion and 
not necessarily when u a — > 0. 


An alternative form for field amplification can be developed through the introduction of the flux 
skin depth, which defines the flux penetration depth during the compression process. The concept of flux 
skin depth is illustrated in figure 1 1 . The flux skin depth is defined by taking the magnetic field at the 
surface as an effective constant value and requiring that the flux contained within a skin depth, , be 
equivalent to the total diffused flux. Conserving the total flux (contained and leakage). 


yielding 




)=MHn ( 


r s s <j>,s 




(26) 


H = j) _ . (2T) 

0 $<p t s ) — ( r a “1“ s <p,a) 

Combining eqs. (22) and (27) gives an expression for A in terms of the flux skin depths in both the armature 
and stator: 

x t> ('•{-'s') 

~{ r a~ s <t>,a) 

For a reasona ble estimation of the compr ession dy namics, it is generally sufficient to use appropri- 
ate approximations or analytical solutions for s^. This approach is known as the skin layer method in which 
it is assumed that the flux loss and penetration velocity of the diffused magnetic field remain limited. That 
is, Sq «{r s - r a ), ds^idt «u a , s^ a « r a , and s « d where d is the thickness of the stator. 
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s e s = Flux Skin Depth in Stator 
s ga = Flux Skin Depth in Armature 


Figure 1 1 . Illustration of flux skin depth concept as applied to flux 
compression reactor geometry. 


4.2 Flux Skin Depth in Planar Geometry 

To determine the time-varying flux skin depth in a semi-infinite planar geometry, the magnetic 
diffusion equation is solved assuming an exponentially increasing magnetic field H/H 0 = x(0 = expt/z c 
where T c is the characteristic time for magnetic compression and is associated with the plasma expansion 
speed. The planar geometry approximation can be invoked whenever flux skin depth is smaller than the 
radius of curvature. 

The magnetic diffusion equation for a semi-infinite planar one-dimensional (distance and time as 
independent variables) geometry with a magnetic field oriented perpendicular to the surface has the form 


1 dH d 2 H A 

+ — t- = 0 , 

D m dt dg 2 

where the following stationary boundary condition applies 

H(0,t) = H 0 e tlTc ■ -oo< t <oo . 

By seeking a solution of the form exp t/z c j{x), a differential equation for/is obtained as 

a 2 / = / 

dg 2 D m r c 


( 29 ) 


( 30 ) 


( 31 ) 
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from which following the general solution is deduced: 


H(g, t) = H 0 e lltc ( V~^° + A 2 e gf ^ ) , ( 32 > 

where 

■ <33) 

The only physically acceptable solution, which satisfies the requirement for a monotonically decreasing 
field in the conductor, is 


H(g,t) = H 0 ex p 


V r c 



(34) 


Again, the flux skin depth is defined such that the diffused flux 0^(0 into the surface can be expressed by 
the surface flux density B SUJ j(t) spread over a depth s ^ : 


0 


(35) 


or 


Hsurf (Oty = J H{g,r)dg . 
0 


(36) 


Thus, evaluation of the integral using the profile defined by eq. (34) yields the desired expression for the 
flux skin depth 


s (p ~ ?0 V ^ni • 


(37) 


4.3 Skin Layer Methodology 

A simple estimation of flux efficiency may be obtained using Sakharov’s skin layer methodol- 
ogy. 31 Although the electrical conductivity will depend on a number of factors which vary during the 
course of compression, consideration is limited to the case of finite and constant electrical conductivity in 
order to avoid the need to solve a system of partial differential equations. The physical motivation for the 
skin layer methodology is based on the fact that the skin layer depth is dominated by the ever-increasing 
field so that the contribution from the previous initial condition is rapidly lost. 
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In the skin layer methodology, the distribution of current in the compression surfaces is approxi- 
mated by an exponential profile: 


s j & a 

Js=J0,s e 9 and Ja=J0,a e 


(38) 


These relationships are exact if H increases exponentially on the boundary; i.e., H/Hq = exp//r t .. It is clear 
that the following differential relationship must then hold true: 


dH _ H 
dt 


(39) 


Re- 


introducing eq. (25), the following is obtained: 


R c = 


H 


dHldt 


( N 

r 2 -r : 
\'s r a J 


dX 

dt 


.-i-I 


(40) 


Assuming the change in flux as a correction factor (dX/dt « 1; A - 1), the following approximation 
is arrived at for the time-dependent variation in T c : 


T c = 


H 


r s Zla 


(41) 


dH, dt 2 r a u a 


This approximation implicitly assumes that the flux loss and penetration velocity of the differential fields 
remain limited (i.e., « r s - r a and ds^/dt « M a ). 

With reference to figures 4 and 6, Farady’s law is applied to the stator and armature to obtain 


_ d( t> + d <t>a 

(42) 

dt dt 


S ’ 

II 

(43) 


dt 


where (f> = (j) Q - (p s - (p a is the trapped flux. Combining eqs. (42) and (43) yields 


d(j) 


c E a dl-fcEa<l\ = — 


(44) 


or 


2KrX+2*r a El=& • 


(45) 
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where it was taken into account the fact that E s and E a are in opposing directions due to Lenz’s law. 
E*=E+(uxB) is the electric field detected by an observer moving with the medium. Therefore, 


77* J0,S 
£\ =— - 


and 


77 * ha 
E n =— - 


(46) 


and eq. (45) may be written as 


2Kr s J ^- + 2nr a ^- = ^- • 
<*s °a dt 


(47) 


Now, Maxwell’s equations require that Vx H = 0 at the compression surfaces such that 


dH = . 

ac ~ J 


(48) 


Using the assumed current density distribution (j = j 0 exp in eq. (48) yields the differential relation- 
ship: 


dH . -{fa 

Tr ,oe 


(49) 


Separating variables and integrating. 


dH = -j 0 l~e V s * d£ 


(50) 


H = -i 


J0 S ((> 


-r'f T” 

? =/ s <t> — 

-0 


Jo s <p 


(51) 


Therefore, eq. (47) may be written in the form 


d(j) _ 2itr s U 2jtr a H 

dt &s s <p,s ®a^<p,a 


By definition, the trapped flux is given by 

0 = flHn(r} - c) 


(52) 


(53) 
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and eq. (52) becomes 


d0 

dt 


r s t r a 


2 0 


^s^<p,s ^a s <p,a 


Hi 


(t-t) 


(54) 


Defining drjdt = ~u a such that dt = -drju a where u a is the velocity of the expanding armature, 
substituting into eq. (54), and separating variables yields 


d(f> 

T 


( 


^ A i& s S(f) s 


+ 


r a 

2dr a 

!^®a s _<p,a j 

Ua[ r s~ r a) 


(55) 


For a reasonable estimation of compression dynamics accounting for magnetic diffusion losses, the 
flux skin depth approximation is applied for a semi-infinite planar geometry to both the stator and arma- 
ture. Then, eq. (55) takes the following form: 


</0 
' <P 


5 | U It 

^s^D-fc J u ai y r .^ ~ r a ^j 


2dr, 


(56) 


Regarding the change in flux as a correction, the time-dependent relationship is utilized for t c as 
defined by eq. (41) such that eq. (56) becomes 


J0 


Id, 


m y s 


+ T n 


T S 


K 


m,ci 


2^2r~dr a 

I (i-Sf 


(57) 


where the magnetic Reynolds number is defined for the armature as R m a = flG a u a r s . 

It is convenient at this point to nondimensionalize the problem in terms of a flux compression 
coefficient A and a dimensionless independent variable <fj: 


A = 0/0o 

$ = r a/ r s 


(podA - d(j) 

hd£ = dr a . 


(58) 

(59) 


25 


Therefore, eq. (56) takes the desired final form 


dX 
' X 


I a tus_ + ^ r r s 


R, 


m,a 


2j2§dt; 


v^M 2 ) 


2 \ 3 / 2 


(60) 


This equation can be numerically integrated to determine the evolution of the flux coefficient A as a func- 
tion of the armature expansion ratio £ 

4,4 Armature Rebound Condition 

With appropriate design, the increasing magnetic pressure due to flux compression will decelerate 
the expanding plasma armature until it is completely stopped and the motion is reversed. If this rebound 
process does not occur, the plasma will impact the reaction chamber wall, most likely with disastrous 
results. Clearly, stator impact is incompatible with repetitive operation and an estimate of the turnaround 
distance is essential. 

At the turnaround point, the assumption is made that the useful energy production per detonation 
has been transformed into magnetic energy. Thus, 


^ \n,t ^m,0 **" ’ 

where W m , is the magnetic field energy at the turnaround point 


(61) 


Wm.<= ■ 

and W m 0 is the original magnetic field energy in the reaction chamber 


(62) 


W m> o = 


-rhI 


(63) 


W k is the kinetic energy production per detonation, and T] is the energy conversion efficiency with which 
the available detonation energy is transformed into a magnetic field energy. V t and V 0 are the flux contain- 
ment volumes at the turnaround point and at the beginning of armature expansion, respectively: 


V t — 7t fy,a) ^ 


(64) 


Vr\ = M 


2 2 ] 
r s -r ( \z , 


(65) 
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and z is the height of the reaction chamber. Furthermore, it is assumed that 


W k = w D p D nrj L z , m 

where w D is the detonation energy per unit mass. Dividing eq. (61) through by W m 0 and evaluating eq. (27) 
at the turnaround point yields the following relationship: 

ff, vlzfl I I 03. . (67) 

Wm -° 

The turnaround radius, r t , is defined as the radial location at which eq. (67) is satisfied. 
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5. MARK I DEVICE 


In the near term, it is anticipated that meaningful technology demonstration, experiments could 
utilize high-explosive detonation charges in order to establish scientific feasibility, since microfusion and 
HEDM technologies are still several years away from fruition. 

The chief difficulty encountered with high-explosive plasma jets is the limited magnetic Reynolds 
numbers that can be achieved. The detonation velocity of high-explosive plasma jets is typically on the 
order of 10 4 m/sec, and this level of performance can only be achieved through special shaped charge 
configurations. 

Experimental investigations of shaped charge explosives as applied to MHD generator applications 
have been conducted in the past. 32 This included high-speed imaging of the plasma jet and measurement of 
its velocity and electricaFconductivity. Irr these tests, It has been observed that a cohesive, hot, rosy- 
colored plasma slug is projected at speeds on the order of 3x 10 4 m/sec. Figure 12 shows the measured 
electrical conductivity in a plasma jet from a 44° cavity charge of unseeded composition^ (C4). Note that 
the measured conductivity in the slug exceeded 10 4 S/m. Inserting representative values for the jet velocity 
and electrical conductivity into the expression for the magnetic Reynolds number, as defined by eq. (19), 
and assuming a reasonable length scale for a practical device, it is found that 

R,„ = Hf,a 0 uok = 4k(10) 7 (10) 4 (10) 4 (10)-' > 1 . < S8 > 



plasma jet of unseeded C4. 
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A practical demonstration device can be developed using a non-nuclear plasma source, based on 
these results. A demonstration device along the lines of the Mark I configuration depicted in figure 1 3 
would be a logical step in this direction. This device is envisioned as a Vz-m-diameter reaction chamber, 
which uses colliding plasma jets from opposing high-explosive charges to produce a radially expanding 
armature. The formation of a radially expanding plasma jet from two colliding jets has been previously 
demonstrated on a small scale. 32 The observed collisional process in these experiments is illustrated in 
figure 14. 


Shielding Magnet 



Explosive 





Figure 13. Mark I configuration for a radial- 
mode explosively driven demon- 
stration device. 


Figure 14. Illustration of plasma jet 

collisional process and radial 
armature formation. 


5.1 Representative Calculations 

As a precursor analysis to the Mark I configuration, calculations based on the skin layer methodol- 
ogy have been performed. The calculations presume an internal generator diameter of 0.5 m and various 
stator materials were investigated. Typical parameters for high-explosive detonations used in the calcula- 
tions are summarized in table 5. 
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Table 5. High-explosive characteristics. 


Density (p D ) 

1,700 kg/m 3 

Specific Energy (w D ) 

5 x 10 6 J/kg 

Velocity (u D ) 

1 xIO 4 m/sec 


The variation in the flux coefficient as a function of the armature expansion ratio is shown in figure 
15. Numerical integration of eq. (60) was performed for magnetic Reynolds numbers of 1 and 10, assum- 
ing stator materials of steel, copper, and a hypothetical HTSC. For the hypothetical HTSC stator material, 
a magnetic diffusivity value of 10 - *’ m^/sec was assumed. With a magnetic Reynolds number of unity, an 
increasing compression efficiency was observed as the material conductivity increases, although the result 
for HTSC is not as encouraging as desired. When the magnetic Reynolds number is increased to 10, a 
significant increase in compression efficiency was noted for all stator materials. 

Calculations for the armature turnaround radius as a function of initial seed field were carried out 
for a copper stator assuming an energy conversion efficiency of 50 percent. The results are presented in 
figure 16 for various magnetic Reynolds numbers. There is a minimum initial seed field, the value of which 
depends on the magnitude of flux diffusion losses, to avoid stator impact. With a good stator material and 
high R m , it is possible to reduce the seed field and minimize system size and weight. 

Apparently, significant flux compression can be accomplished at marginal magnetic Reynolds num- 
bers when the stator material is sufficiently resistant to magnetic diffusion. Thus, it appears that a success- 
ful demonstration device could be developed, which utilizes a non-nuclear high-explosive plasma jet source. 
The preceding analysis is highly idealized, of course, and real hydrodynamic effects will lead to more 
pessimistic predictions with respect to the minimum required magnetic Reynolds number; however, the 
simplified analysis results appear encouraging. 
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6. MAGNETIC DIFFUSION LABORATORY EXPERIMENTS 


6.1 Analytical Model 

As part of this research, simple laboratory experiments were designed to investigate the magnetic 
diffusivity characteristics of HTSC samples when exposed to strong-pulsed magnetic fields. The basic 
configuration in these experiments, as illustrated in figure 17, is a hollow tube of test material surrounded 
by a solenoid, which can be used to create a pulsed external magnetic field on demand. Measurement of the 
time-varying magnetic field inside the tube, once the solenoid is pulsed, yields quantitative information 
with respect to the magnetic diffusivity properties of the material. In sections 6.1 and 6.2, an analytical 
treatment of the magnetic diffusion process through a conducting hollow cylinder is described and the 
preliminary experimental results discussed. 



Figure 17. Illustration of experimental configuration 
for examining pulsed magnetic diffusion 
through a hollow conducting cylinder. 


The magnetic diffusion through a conducting hollow cylinder is governed by a magnetic diffusion 
equation derived from Maxwell’s electromagnetic laws. In a cylindrical coordinate system, the governing 
diffusion equation has the form 33 


d 2 H 1 dH ? 1 dH, 
dr 2 r dr D m dt 


( 69 ) 


33 


and the following initial/boundary conditions apply: 33 



H z (b,t) = H 0 

(70) 


H z (r, 0) = 0 

(71) 

dH z 

_2D /n ra^i . 

(72) 

dt 

a a [dr \ r=a 



The solution to this system of equations is known from mathematical physics using the Laplace 
transform methodology. The field H z (t) = //_(<:/,/), which is built up in the cavity of the conductor, can be 
given for fl R = 1 in the form 33 


H z (a,t) = H 0 -4H 0 ^e- 

n = 1 


J 2 (aa n )J 0 (bet A ) 


(aa n f Jl(ba n )-Jl (««„)] 


where the a n are roots of 

jQ(ba)Y 2 (aa ) - YQ(ba)J 2 {aa) = 0 

and J n and Y n are Bessel functions of the first and second kind of order n, respectively. 
After a long enough time, such that 


(73) 


(74) 


t> j ’ 

D m (X\ 

it is sufficient to retain only the first (n = 1) term. The field growth in the conductor is then 
where 

1 _ b > 2 

T ‘ D m al 4 nD m 8 


(75) 


(76) 


(77) 
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The parameter g, defined by 


4k 

8 ~ (Z?aj) 2 ’ 

is a geometrical form factor which has been tabulated as a function of the ratio a/b in table 6. 34 


( 78 ) 


Table 6. Tabulation of geometric form factor 
as a function of the tube radius ratio. 


a/b 

9 

a/b 

9 

0.00 

2.173 

0.65 

1.587 

0.05 

2.173 

0.70 

1.431 

0.10 

2.172 

0.75 

1.253 

0.15 

2.170 

0.80 

1.052 

0.20 

2.164 

0.85 

0.8267 

0.25 

2.153 

0.90 

0.5766 

0.30 

2.133 

0.95 

0.3011 

0.35 

2.102 

0.96 

0.2430 

0.40 

2.059 

0.97 

0.1838 

0.45 

2.001 

0.98 

0.1236 

0.50 

1.926 

0.99 

0.0623 

0.55 

0,60 

1.833 

1.720 

1.00 



In our small-scale experiments, the tube size was as follows: length = 75 mm, wall thickness 
= 7. 14375 mm, and outside diameter = 30. 1 626 mm. This yields an inner radius of a = 0.079375 cm and an 
outer radius of b = 1 .508125 cm. Thus, g = 1 .8771 (from table 6) and the characteristic diffusion time may 
be estimated according to eq. (77). The characteristic diffusion times, assuming a range of conductivity 
values for the HTSC material, are summarized in figure 1 8. Actual diffusion characteristics of the HTSC 
material will depend on the behavior of the superconductor under strong pulsed fields. Based on the 
experimental efforts of previous researchers, D m is anticipated to fall in the range of 10 -5 to 10 -3 m 2 /sec. 
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Magnetic Diffusion Time Constant (sec) 



Magnetic Diffusivity (m 2 /sec) 


Figure 18. Characteristic magnetic diffusion time constant for magnetic 
field penetration through conducting hollow cylinders of 
various materials with g = 1.8771. 
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6.2 Laboratory Experiments 


The objective of the experiments was to measure the magnetic diffusion rate of Bi 2 Sr 2 CaCu 2 O x 
(BSCCO) superconducting material in comparison with conventional metals. 

The experimental apparatus consists of a cylindrical tube surrounded by an excitation coil. A Hall 
probe is placed along the centerline of the tube to measure the instantaneous magnetic field, and a Rogowski 
coil is looped around the tube and coil as a means of measuring the induced current. The entire test section 
is submerged in liquid nitrogen for testing. Different materials (i.e., copper, aluminum, and steel) were 
tested in order to compare their magnetic diffusion rates to that of the superconductor. 

During testing, it can be expected that a time delay will occur between the external excitation field 
and internal field shielded by the tube. From this experimentally measured time delay, the magnetic diffu- 
sion characteristics of the material can be deduced. 

Some typical results from these experiments are provided in figure 19 for aluminum and BSCCO 
materials. These results indicate a time delay of =45 msec for BSCCO which may be compared with a time 
delay of 4 msec for aluminum. Although preliminary in nature, these data indicate a significant improve- 
ment in penetration delay time and are sufficiently encouraging to warrant further investigation. 


Magnetic Field Versus Time for BSCCO Superconductor and Aluminum Tubes 



Figure 19. Magnetic penetration time delay characteristics 
for BSCCO and aluminum tubes. 
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7. CONCLUSIONS 


Magnetic flux compression reactors are technically feasible for spacecraft propulsion and power, 
based on our preliminary assessment. It appears that high-power density systems can be designed which 
can utilize a detonation plasma armature in conjunction with a highly conductive stationary stator. In this 
way, one could develop an intermittently fired reactor using microfusion or advanced HEDM targets for 
the direct production of electrical power and thrust. 

A Mark I demonstration radial-mode reactor driven by chemical high explosive is also proposed. 
Analyses of the conductivity and velocity produced by shaped charge plasma jets indicate that one can 
achieve R m > 1 using non-nuclear sources. A skin layer analysis of the concept indicates that effective flux 
compression can be achieved in a relatively moderate scale device. In addition, the potential for improved 
field compression using type II HTSC stator materials was evaluated, but the performance gains were 
marginal based on the assumed conductivity value. Small-scale laboratory experiments did show that BSCCO 
has a significantly lower magnetic diffusivity in comparison to conventional metals, but further research is 
needed to clarify the actual system-level benefits. The major technical issues associated with the use of a 
plasma armature appear to be achieving good diamagnetic properties and controlling Rayleigh-Taylor 
surface instabilities. 
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